4200 OPTICS LETTERS / Vol. 37, No. 20 / October 15, 2012 


Infrared dielectric properties of low-stress silicon nitride 
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Silicon nitride thin films play an important role in the realization of sensors, filters, and high-performance circuits. 
Estimates of the dielectric function in the far- and mid-IR regime are derived from the observed transmittance spec- 
tra for a commonly employed low-stress silicon nitride formulation. The experimental, modeling, and numerical 
methods used to extract the dielectric parameters with an accuracy of approximately 4% are presented. © 2012 
Optical Society of America 
OCIS codes: 310.3840, 310.6188, 310.6860. 


The physical properties of silicon nitride thin films, 
namely low tensile stress, low thermal/electrical conduc- 
tance, and its overall compatibility with other common 
materials, have facilitated its use in the microfabrication 
of structures requiring mechanical support, thermal iso- 
lation, and low-loss microwave signal propagation (e.g., 
[1-4]). Silicon nitride films are amorphous, highly absorb- 
ing in the mid-IR [5], and their general properties are 
functions of composition [6,7]. Here the optical proper- 
ties are studied in detail for a membrane with parameters 
commonly employed in microfabrication. 

The silicon nitride optical test films were prepared by a 
low-pressure chemical-vapor-deposition (LP-CVD) pro- 
cess optimized for low tensile stress and refractive index 
[8]. The 5 ! 1 SiH 2 Cl 2 /NH 3 gas ratio employed results in a 
tensile stress <100 MPa and optical index greater than 
~2 [9]. The test structure is shown schematically in Fig. 1 
(inset). Double-side-polished silicon (75 mm diameter, 
500 pm thick) wafers [ 10 ] were used as a mechanically 
robust handling structure for the SiN,.. membranes. A 
150 nm thermal oxide was grown on the silicon wafers 
by wet oxidation at 950°C for 31 min. This layer was sub- 
sequently used as an etch stop to protect the nitride dur- 
ing definition of the silicon handling wafer geometry. A 
low-stress SiN,. layer was then deposited by LP-CVD [e.g., 
deposition parameters for 2 pm film are 835°C for 9.7 h 
with pressure 33 Pa and 12 seem NH3, 59 seem SiH 2 Cl 2 
(SCCM denotes cubic centimeters per minute at standard 
temperature and pressure)]. The wafers were then pat- 
terned with a resist mask and SiN r /Si() 2 windows 
formed by deep reactive ion etching, which removed 
all the silicon under the window area. The residual 
thermal oxide was removed with hydrogen fluoride va- 
por etch, leaving a set of uniform SiN,, membranes each 
with a 10 mm diameter aperture individually suspended 
from the silicon handling frame. 

The optical tests were performed on SiN r samples hav- 
ing membrane thicknesses of 0.5 and 2.3 pm with a un- 
certainty of 3%. Fabry-Perot resonators were made by 
stacking multiple samples with silicon standoff frames 
between adjacent samples to explore the long- 
wavelength response of the material in greater detail. 
The silicon standoffs allowed a vent path for evacuation 
of air between the nitride membranes. All optical 


measurements were performed in vacuum with a resi- 
dual pressure less than 100 Pa. 

The samples were characterized with a Bruker 125 high- 
resolution Fourier transform spectrometer (FTS) and 
were measured in transmission at the focal plane of an 
/ /6 beam. A number of different sources, beam splitters, 
and detector configurations were used in combination to 
provide measurements over the reported spectral range. 
The single-layer SiN,, sample transmission was measured 
over an extended range from 15 to 10,000 cm -1 . The mer- 
cury lamp and a multilayer Mylar beam splitter were used 
to access frequencies below 600 cm -1 . Additional mid-IR 
spectral data up to 2400 cm -1 were acquired using a cera- 
mic glow bar source, Ge-coated KBr beam splitter, and 
room-temperature deuterated tri-glycine sulfate detector. 
The remaining near-IR data up to 10, 000 cm -1 were taken 
with a W filament source, Si on CaF 2 beam splitter, and a 
liquid-nitrogen-cooled InSb detector (Fig. 1). Far-IR data 
between 15 and 95 cm -1 were taken using a mercury arc 
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Fig. 1. (Color online) Room-temperature transmission of a si- 
licon nitride sample 0.5 pm thick: measured (grey), model 
(black dotted), and residual (red). The shaded band’s width de- 
limits the estimated 3<r measurement uncertainty. A 30 GHz 
(1 cm -1 ) resolution is employed for the measurement. The in- 
sert depicts the geometry of the SiN, membrane and microma- 
chined silicon frame. 
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Fig. 2. Measured (solid grey) and model (black dotted) trans- 
mission for a three-layer stack of silicon nitride samples 2.3 pm 
in thickness with 998 pm intermembrane delays that comple- 
ments the data shown in Fig. 1. The sample response in the 
far-IR was acquired with a resolution of 3 GHz (0.1 cm -1 ). 


lamp source and a liquid-helium-cooled 4.2 K bolometer. 
Mylar beam splitters of 50, 75, and 125 pm thicknesses 
and a multilayer Mylar beam splitter were used during se- 
parate scans (Fig. 2) . The resultant transmission data were 
merged into a single spectra using a signal-to-noise weight- 
ing for subsequent parameter extraction. 

The dielectric response is represented as a function of 
frequency, w, by the classical Maxwell-Helmholtz-Drude 
dispersion model [11], 


£,.(&>) = £«, + 


cw 2 . - co 2 - iwT'j(w) 


( 1 ) 


equation to the laboratory FTS data. Specifically, a se- 
quential quadratic programming method with computa- 
tion of the Jacobian and Hessian matrices [14,15] was 
implemented. The merit function, y 2 , was used in a con- 
strained minimization over frequency as follows: 

J2 h) - TVrsJ 2 , (3) 

k= 1 

where N is the number of data points, T the modeled 
transmittance, 7Vrs the measured transmittance data, 
and h the measured sample thickness. We are guided 
by the Kramers-Kronig relations in defining constraints 
for a passive material: Cy > £y +1 1, e'j > 0 and £ r (0) = iq 
[16], For accurate parameter determination, the sample 
should have uniform thickness, be adequately transpar- 
ent to achieve high signal to noise, and have diffuse scat- 
tering as a subdominate process. The method requires an 
a posteriori numerical verification for Kramers-Kronig 
consistency. In the example presented here, a numerical 
Hilbert transform [17] of e' r (oo) reproduces e" (o>) to with- 
in 2% (Fig. 3). An alternative method employing reflectiv- 
ity and phase allows a priori Kramers-Kronig consistent 
results [18]. However, given the details of the thin-film 
samples and available instrumentation, this approach 
was not implemented. 

Figure 1 illustrates the measured and modeled results 
obtained from the analysis of a 0.5 pm thick sample. The 
peak residual in the transmittance is less than 3%, and the 
3o — 0.023 uncertainty band indicated corresponds to 
the 99.7% confidence level. The standard deviation 
adopted for the measured data, a, was estimated assum- 
ing the errors as a function of frequency were uniform 
and had a reduced y 2 equal to unity. An additional 


where M is the number of oscillators and e r — e' r + ie". is 
a complex function of (5 M + 2) degrees of freedom, 
which are as follows: the contribution to the relative per- 
mittivity e„ = £,v + i of higher lying transitions, the differ- 
ence in relative complex dielectric constant between 
adjacent oscillators A ej — £y - £y +1 , which serves as a 
measure of the oscillator strength, the oscillator reso- 
nance frequency oo T ., and the effective Lorentzian damp- 
ing coefficient T'-, for j — 1 M. The following 

functional form is used to specify the damping: 


P(cy) = Pexp 


a>r - OOt 

^7 


,2 X 2 n 


( 2 ) 


where aj allows interpolation between Lorentzian 
(ay — 0) and Gaussian wings («y > 0) similar to the ap- 
proach in [12]. The form indicated above enables a more 
accurate representation of relatively strong oscillator 
features. 

The impedance contrast between free space and the 
thin-film sample forms a Fabiy-Perot resonator. The 
observed transmission can be modeled [13] as a function 
of the dielectric response [Eq. (T)], thickness, and wave- 
number. The dielectric parameters were solved by means 
of a nonlinear least-squares fit of the transmission 
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Fig. 3. (Color online) Real and imaginary parts (solid red 
curves) of the dielectric function as extracted from the data 
shown in Fig. 1. The line thickness is indicative of the propa- 
gated ~4% error band. The numerical Hilbert transform of 
the modeled e" (o>) is indicated in the upper panel (dashed blue 
line) to facilitate comparison with e' r (co). The filled symbols 
indicate the parameters derived from the data presented in 
Fig. 2. 
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Table 1. Fit Parameter Summary 


3 h] 

<■ H 

e 'i H 

co Tj /2n [THz] 

Vj/2ji [THz] 

a j H 

1 

7.582 

0 

13.913 

5.810 

0.0001 

2 

6.754 

0.3759 

15.053 

6.436 

0.3427 

3 

6.601 

0.0041 

24.521 

2.751 

0.0006 

4 

5.430 

0.1179 

26.440 

3.482 

0.0002 

5 

4.601 

0.2073 

31.724 

5.948 

0.0080 

6 

4.562 

0.0124 





uncertainty in the FTS normalization influences the 
dielectric response function at the 1% level. In addition 
to the channel spectra, the observed spectrum shows 
two predominant features at 12 and 25 THz. Simulations 
with M — 2 oscillators lead to a peak residual on trans- 
mission of 5% and do not enable recovery of the reso- 
nance at 25 THz. Using five oscillators satisfactorily 
recovers the observed transmittance and reduces the 
peak residual by a factor of 4.4. When the resonator’s 
quality factor, Q eS — a T /Tj, is greater than 5, the data 
were not reproducible by either a pure Lorentzian oscil- 
lator or Eq. (4.6) in [12]. In these regions, the peak trans- 
mission residuals were decreased by a factor ~2 through 
the use of Eq. (2). 

In Fig. 3 the values of the real and imaginary compo- 
nents of the dielectric function are illustrated as a func- 
tion of frequency. The uncertainty in e r was propagated 
and computed as described in [19]. Table 1 contains a 
summary of the best fit parameters for five oscillators, 
which can be used to reproduce the data shown in 
Fig. 3. 

To characterize the long-wavelength portion of the di- 
electric function, Fabry-Perot resonators were realized 
from one-, two-, and three-layer samples. Representative 
data for the three-layer resonator stack are presented in 
Fig. 2. A multilayer transfer matrix analysis [13] is used to 
extract the dielectric function using the measured SiN,. 
(2.3 pm) and silicon spacer (998 pm) thicknesses. The cir- 
cular symbols at 1.5 and 2.5 THz indicated in Fig. 3 were 
computed from a composite analysis of the three Fabry- 
Perot measurement sets. The horizontal range indicates 
the data used in each fit. The best estimates are e r ~ 7.6 + 
i0.08 over the range of 2-3 THz and e r ~ 7.6 + 70. 04 
over 0.4-2 THz. The real component of the static dielec- 
tric function derived from the data is in agreement with 
prior reported parameters for this stoichiometry [4], As 
shown in Fig. 3, the measurements are internally consis- 
tent and represent roughly a factor-of-three reduction in 
uncertainty relative to prior IR SiN r; measurements 
identified by the authors [5-7]. The dielectric parameters 
reported here are representative of low-stress SiN r 


membranes encountered in our fabrication and test 
efforts. 
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